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The transresistance and the Hall transresistance of dirty two-dimensional bi-layer systems in the 
presence of tunneling bridges (pinholes) are studied theoretically. We find, at weak magnetic field, 
a non-zero Hall transresistance. In a geometry where the pinholes are concentrated in the middle of 
the sample, a quantum process gives the dominant contribution to both the ordinary transresistance 
and the Hall transresistance. Arising from the interplay between Coulomb repulsion, disorder and 
tunneling, the quantum contribution increases in a singular way as the temperature decreases. 
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I. INTRODUCTION 

The progress in micro-structure technology of semi- 
conductors has made it possible to fabricate a pair of 
parallel two dimensional (2D) electronic layers that are 
spatially close to each other. Experimental and the- 
oretical efforts p|-|l5[ are being performed to understand 
the behavior of these systems. The physics of bi-layer 
system is interesting in its own right , but it also serve 
as a tool to test the internal layer properties. 

In a typical experimental set up, a current / is sent 
through one of the layers (layer 1) and a voltage drop, 
Vt, is measured by separate contacts on the other one 
(layer 2) (see Fig. [j]) . The ratio between Vt and / defines 
the transresistance, Rt- In a similar way we can define 
a Hall transresistance , when a magnetic field is ap- 
plied. [For precise definitions see Eqs. (|35|) and (|36| ) in 
Sec. [0.] The behavior of the transresistances is rather 
rich due to the possibility to control the layer areas, the 
electron density in the 2D layers, their mobility, the in- 
terlayer tunneling rate, and external parameters like the 
temperature and a magnetic field. 

In this work we concentrate on the corner of the param- 
eter space where the mobility is relatively low (disorder is 
large), a tunneling between the layers through local pin- 
holes (or bridges) is possible and a weak magnetic field 
can be turned on. This situation occurs when the aver- 
age distance between the layers, d, is relatively small so 
we expect the Coulomb forces to be dominant |ll|. [We 
consider here only the cases of a weak or zero magnetic 
field, with weak interlaycr tunneling and weak interlayer 
interaction; generalizations to strong magnetic fields and 
nonperturbativc tunneling and interactions are subjects 
for future studies.] 

In the absence of tunneling, a transresistance arises 
by frictional (drag) forces due to the Coulomb repulsion 
between electrons in the two layers. This mechanism in- 
volves thermal density fluctuations around a mean value 
of the electron density and therefore vanishes at low tem- 
peratures. In the remainder of the article, this effect will 
be referred to as the classical drag mechanism. 

In the presence of tunneling through pinholes there are 



additional physical processes that lead to a finite tran- 
sresistance. The first has to do with the fact that, in 
the presence of pinholes, current can flow from layer 1 
to layer 2 through a pinhole close to the source lead and 
flow back through another pinhole close to the drain lead. 
This purely classical effect leads to a net current flow in 
layer 2, to a potential drop and finally to finite transresis- 
tance, which we refer to as the leakage contribution. The 
potential due to this mechanism depends on the distribu- 
tion of the pinholes in the sample and, through the weak 
temperature dependence of the pinholes' conductance, on 
the temperature. 

In addition, there is a more subtle mechanism due to 
the interplay between Coulomb repulsion, tunneling and 
disorder |l5[ |. We shall see that the main part of this ef- 
fect arises from frequencies that are larger than the tem- 
perature [see discussion after Eq. (^3|)]. In that sense, 
the effect involves virtual processes of quantum origin, 
and will be called the quantum mechanism. The quan- 
tum effect is a generalization of the intralayer interaction 
correction in a single layer disorder system to the 
geometry of bi-layer systems. 

Unlike the classical drag mechanism, which vanishes 
at low temperature, and the leakage contribution, which 
depends weakly on the temperature, due to the electrons 
diffusive motion in layer 1 and layer 2 the quantum mech- 
anism increases in a singular way when the temperature 
decreases. 

Without pinholes the Hall transresistance vanishes 
rapidly (cx T 4 ) at low temperatures |ll],[l8|]. (We do not 
discuss the case where there are strong correlations be- 
tween the layers ITS].) We will see below, however, that in 
the presence of tunneling, there are nonvanishing leakage 
and quantum contributions to the Hall transresistance. 
Thus, a measurement of the Hall transresistance is a di- 
rect measurement of the leakage and quantum effects. 

In most geometries, it might be difficult experimen- 
tally to separate the quantum and the leakage contribu- 
tions because the temperature dependence of the tran- 
sresistance measurements has to be very accurate [ p"5| . 
However, we argue below that in a geometry where the 
pinhole distribution is deliberately concentrated near the 
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middle of the sample, at low enough temperature the 
quantum mechanism is larger than both the classical 
drag mechanism and the leakage contribution. While the 
standard drag measurement in the absence of tunneling 
gives information on thermal fluctuations and interlayer 
interactions of the system, in the presence of tunneling 
through pinholes the transresistance measurements can 
provide interesting information about quantum processes 
that involve an interplay between disorder and interac- 
tion. 




distribution is given in Sec. [v| below, for several cases of 
interest. 

In case where L2 3> W2 and the measuring probes for 
the longitudinal transresistance R t are very far from the 
tunneling places, the potential in layer 2, (x, y), is 
practically independent of y near the points of measure- 
ment. In that case the average over y in the definition of 
Rt [see Eq. (p^) ] is not needed. However, in this geometry 
the Hall transvoltage must be measured close to the tun- 
neling point and is sensitive to the distance from it along 
the x axis. [In fact, as shown in Sec. V A, the Hall trans- 
voltage falls off exponentially with the distance from the 
tunneling points.] To avoid this factor, it is preferable 
to measure the Hall transvoltage in a "cross" geometry 
where L2 ~ W\ <C W2. In that case the Hall transre- 
sistance is measured far from the tunneling points, the 
potential V^ 2 \x, y) depends weakly on x, and the aver- 
age over x in Eq. ( pq ) is not needed. 

We shall see, below, that Rt and i?f may be written 
as the x and y components of a vector R tl which has the 
form 
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FIG. [l] Geometry for a drag experiment. The lightly shaded 
areas denote the overlapping regions of the 2D electron gases. 
In a typical transresistance measurement a current I is sent 
through the 2D layer 1 and a (trans)potential Vt is developed 
in the 2D layer 2. In the absence of tunneling between the 
layers, the Hall transpotential V t H is zero at a weak magnetic 
field H perpendicular to the layers. In the text we discuss 
the influence of tunneling through local pinholes (or bridges) 
between the layers on the transresistance and the Hall tran- 
sresistance [defined in Eqs. (^) and (fjq)], in the presence of 
Coulomb repulsion. The dark shaded regions denote the area 
where pinholes exist and tunneling between the layers is pos- 
sible. We consider cases where the tunneling region is much 
smaller then the overlapping region and where the two regions 
are equal. 

For simplicity we quote here results for 2D layers that 
have identical properties; i.e., they have the same sheet 
resistance R\j, diffusion constant D, Fermi energy Ep, 
Fermi momentum kp, mean free path I and mean free 
time t, total density of states (including spin) v and in- 
verse Thomas- Fermi screening length n — 2ire 2 u. 

Since the measured transvoltages V t and V t H depend 
on the precise location of the voltage contacts, we use 
here an average over the voltages along the appropriate 
edges in the definition of the transresistances; e.g., for 
rectangular geometry we integrate the potential along the 
boundary and divide by its length. [For precise definition 
of the transresistances see Sec. [II C| .] We assume that 
the length L\ of layer 1 is large compared to its width 
Wi, so that the current density is uniform through the 
interaction region. A discussion of the actual voltage 
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where I is the total current flowing in layer 1, and W2 
is the width of layer 2 (see a gener alization for nonrect- 
angular configurations in Sec. HI C). The components of 
the vector P, arising from the leakage contribution, are 
given by the product of resistivity tensor of layer 2 and 
the dipolc moment of the tunneling current distribution. 
The vector F is due to the momentum transferred from 
layer 1 to layer 2, it includes both the classical drag con- 
tribution and the quantum effect in the presence of tun- 
neling between the layers. 

If the tunneling between the layers occurs uniformly 
all over the sample then the dipole moment P is large 
and the quantum mechanism is a small correction to the 
leakage contribution to the transresistances. However, 
when the pinholes are concentrated in the middle of the 
samples the dipole moment is small and the quantum 
contribution is dominant. 

When layer 1 and layer 2 have rectangular shapes of 
sizes L\ x W\ and L2 x W2 respectively and the pinholes 
distribution is concentrated in a rectangle of dimensions 
a x b near the middle of the sample we find: 
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where po oc T 2 is the bulk drag coefficient [see the pre- 
cise expression for it in Eq. (P7|)1, in = i?ne 2 /2tt 2 H, 
Si a t — L x W is the area of the layers' overlapping re- 
gion, L = min{Li,L 2 } and W = min{Wi, W2}, Rh is 
the Hall resistivity of a single layer, and R± is the total 
resistance between the layers. The term proportional to 
a 2 (or a 2 +b 2 for the Hall transresistance) is from the leak- 
age contribution and the term proportional to L 2 ^ — D/T 
is the quantum contribution. The coefficients at and ol^ 
are numbers of order unity and the corresponding terms 
are due to the zero bias anomaly correction to R±, and 
due to intralayer interaction and weak localization cor- 
rections to the conductivity within each layer j2(J (see 
also Sec, gvg ). 

Expressions and (|J) are valid for 
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where L m ; n = min-fL,]/^}, H is the external magnetic 
field perpendicular to the layers, and p is the sample 
mobility. 

If L2, W2 and W\ are all comparable to each other, 
then the quantum contribution simply flattens out and 
becomes temperature independent for T < D/L 2 nhl . 
On the other hand if L 2 is much larger than W± and 
W2, then there could be an intermediate temperature 
D/L\ <C T <C D/L^^ where the system is quasi one- 
dimensional, and the temperature dependence of the 
quantum contribution to Rt may be even more singu- 
lar then in Eq. (0). If R± is not very large so that the 
energy scale ( 1 /r) exp [~n Rj_/R^j] may be larger than 
-D/Lj^jn, then effects which are non linear in RJ 1 may 
need to be taken into account at low temperatures. If 
H exceeds 1/p then effects nonlinear in magnetic field 
become important. 

We have also assumed through out the paper that the 
current in layer 1 is so low that the cut off of the quantum 
process is determined by the temperature and not by 
the voltage difference between layer 1 and layer 2. This 
assumption should hold as long as JR^jLt, J Rub «T/e 
where J is the current density in layer 1. 

Let us examine now expressions (||) and (^) for the 
transresistances. If we further assume that a,b <C Lt <C 
L m in then the leakage contribution is suppressed and the 
quantum contribution is larger than the leakage contribu- 
tion. On the other hand, the classical drag contribution 
could be larger than the quantum contribution at finite 
temperatures, even though pp vanishes as T 2 , because 
the classical drag is effective over the entire area of the 
overlap region, as reflected in the prefactor Si nt . In order 
to minimize the classical drag contribution, one should 
choose the dimension of the sample as small as possible, 
consistent with the requirement that L m ; n remain larger 
than Lt at the lowest achievable temperatures. 

By contrast, there is no contribution from the classical 
Coulomb drag to the Hall transvoltage. This happens be- 
cause in the absence of tunneling no current is flowing in 



layer 2, there is no Lorentz force that should be compen- 
sated, and as a result no Hall transvoltage is developed 
at low temperatures Q . 

In case when a, b are comparable to the layer size, i.e., 
when the pinhole distribution is uniform, the quantum 
contribution is a small correction to the leakage contri- 
bution, similar to the small interaction corrections in sin- 
gle layer substances ^(J. In that case the temperature 
dependence is mainly determined by the intralayer inter- 
action, and weak localization corrections to R\j and the 
zero bias anomaly correction to R± . 

Examining Eq. (||) we see that at low temperatures the 
contribution from the classical drag vanishes and the con- 
tribution from the leakage term oc a 2 is a temperature- 
independent constant. The quantum contribution in- 
creases as 1/T. Thus with the right choice of parameters 
and at low enough temperatures the quantum contribu- 
tion is dominant. We note that while the usual drag 
mechanism leads to a potential drop in layer 2 that is 
opposite to the voltage drop in layer 1, both the leak- 
age and the quantum mechanisms give rise to a potential 
drop in the same direction as in layer I. Therefore pn has 
a sign which is opposite to leakage and the quantum con- 
tributions. Thus we expect pl| ] that at a temperature 
T* the transresistance will change signs. The Hall tran- 
sresistance has no contributions from the classical drag 
mechanism JO], hence, it gives a direct measurement of 
the leakage and the quantum contribution. 

For a GaAs /AlGaAs rectangular sample of dimensions 
20pm x 5pm, mobility fi = 5x I0 4 cm 2 /Vs, electron den- 
sity n = 4x 10 10 cm~ 2 , i? ± = 20fcSl and a = OApm, b = 
0.1/im, yielding R\j = 3fci7, we find (with nd = 3) that 
the total contributions of the classical, leakage and quan- 
tum mechanisms (neglecting zero bias anomaly, weak lo- 
calization and intralayer interaction corrections) are: 
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with T measured in Kelvin, for T > 2mK. At T* ^ 0.3K 
the transresistance is zero. This means that the leakage 
and quantum contributions win over the classical drag 
and at T < 0.2K the quantum contribution is larger 
than the leakage one. At T = 2mK the system become 
quasi-I-D and the behavior of the quantum corrections 
is even more singular cx I /T 3 / 2 [See Eq. ( 49 )] , eventually 
the at T = QAmK the quantum contribution becomes 
temperature independent. 

The Hall transresistance in the "cross" geometry, i.e., 
when the dimensions of layer I are 20pm x 5pm and of 
layer 2 are 5pm x 20pm with the same 2D electron gas 
parameters as above, is found to be 



R? (mO) S (160 + ^)H, 
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with H measured in Tesla, for T > 2mK, H < 0.2T. 
The last condition insures that we are in the linear regime 
with respect to the magnetic field. Since in the absence of 



3 



tunneling the Hall transresistance is zero, this finite Hall 
transresistance is a direct measurement of the leakage 
and quantum corrections. 

We note that the tunneling region in the middle of the 
sample does not have to be a square. If it has the shape of 
a slit geometry, e.g., a = 0.01pm and b = lpm then the 
leakage term in Rt is even smaller. One should have in 
mind, though, that in this case the slit has to be aligned 
very precisely perpendicular to the current in layer 1, in 
order to keep the leakage term small. In the slit geometry 
it is easier to make R± larger. 

The structure of the remainder of the paper is as fol- 
lows: in Sec. [n]we discuss the formulation of the problem 
in terms of the resistivity tensor p^(r, f), the appro- 
priate boundary conditions and the continuity equation. 
This leads to a set of integro-differential equations (0), 
(||) and (||). Their solution, combined with the appro- 
priate generalization of Gauss's law, gives the transre- 
sistances in terms of the conductivity tensor a"J 3 (r,r > ), 
inverse of the resistivity tensor, which can be determined 
using a Kubo formalism. In Sec. |^ a microscopical anal- 
ysis of different parts of CT,™' 3 (r, f) is performed using 
the linear response formalism (which we generalized to 
include the tunneling through local pinholes) for dirty 
materials. Later we discuss in some m ore d etails the po- 
tential distribution in layer 2. In Sec. VA we solve the 



integro-differential equations in the case where the tun- 
neling occurs only in the middle. We find that the current 
flow lines in la yer 2 are similar to the field lines of a dipole 
in 2D. In Sec. VB we discuss the case when the pinhole 
distribution is uniform. We solve Eq. (Q), (||) and (||) 
perturbatively and find expression for the potential in 
layer 2. Finally, after a concluding section, appendixes 
with details of calculations are presented. 
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where a!-™ are the conductivities of layer a in the absence 

of the other layer and the vector f^°^ , describing momen- 
tum transfer from layer /? ^ a to layer a, is defined be- 
low. To complete the set of the differential equations (Q) , 
(||) we impose the following boundary conditions, (i) No 
current can flow perpendicular to lateral edges of layer 1 
or to the boundaries of layer 2; (ii) the current enters and 
leaves layer 1 with a uniform current density J. When 
layer 1 has the shape of a rectangular of length L\ and 
width W\ , the boundary conditions on Jf (r) are 
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(8) 



j£\±Lx/2,y) = J, A l \x,±Wx/2) = 0, 



and 



(2) 



dS 2 



0, 



(9a) 



(9b) 



where 5*2 is the region of layer 2, dS% its boundary and 
n is a vector normal to dS2- Notice that, since charge 
can not be accumulated in layer 2, (9b|) is possible only if 
JJ S2 ePrV ■ jt 2 ) = 0. A solution of (?])-(|) gives the tran- 

sresistances in terms of f( a \ [Different boundary condi- 
tions reflecting different experimental configurations can 
be analyzed in a way similar to the one discussed below, 
they will change the results for the measured transvolt- 
ages.] 

The components of the vector are given by 
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where 



II. MACROSCOPIC EQUATIONS 

To be specific we will analyze a system with the ge- 
ometry depicted in Fig. [j]. We use a notation where the 
indices i,j, k run over the directions x, y and a, (3 = 1,2 
are layer indices. Summation over repeated indices is un- 
derstood. The local current (per unit length) in layer a 
and direction i, J", is related to the potential difference 
between the layers by the continuity equation: 



ViJf (r) = {-l) a J T if) = gH^A^V^r), 



(7) 



where Jt (r) is the tunneling current density between the 
layers, the matrix A a/3 is —1 if a = (3 and 1 if a ^ j3, 
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and t l is the transmission amplitude of the Zth pinhole lo- 
cated at R l . In addition, the current at point r is related 
to the electric fields via a generalized Ohm's law that 
includes the momentum transfer from the other layer: 



Pkjirf) = pn'5 af }5(r- f) - p^Ur.r 1 ). 



(11) 



The symbol p^ji^f) is a matrix in the variables r and 
f , and in the layer and Cartesian indices, which is the 
matrix-inverse of the conductivity tensor af- (f, r*) that 



can be found from a microscopic treatment (see Sec. IV)] 



The tensor pj^ S(r — f f ) = p^ff, r 1 ) is the resistivity ten- 
sor of layer a in the absence of the other layer (inverse 
to <j( Q )), and S/3 denotes the region of layer (3. The elec- 
tric field (r) describes the field formed in layer a due 
to processes that involve the other layer. For the case 
of weak coupling between the layers, which we consider 
here, only the elements of p which are off-diagonal in the 
layer index need be included, as the diagonal elements 
are negligible. 

There are two essential contributions to p: one is due 
to frictional forces (the standard classical drag) and the 
other is related to the quantum mechanism mentioned 
earlier. [The leakage contribution is captured by the con- 
tinuity equation (^|).] The drag mechanism does not in- 
volve any tunneling through local pinholes between the 
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layers. In this mechanism the electrons in layer 2 interact 
with thermal fluctuations of the electrons current density 
in layer 1 via Coulomb forces. They rectify them and in 
this way are dragged in the direction of the current in 
layer 1. This process leads to accumulation of charges 
at the edges of layer 2. As in the case of a standard 
Hall effect, the charge accumulated at the edges develops 
an electric field that opposes and cancels the force on 
the electrons due to drag. Therefore the voltage drop in 
layer 2 is opposite to the direction of the current in layer 
1, i.e., the drag contribution, pp, to the transresistance, 
Rt, has a negative sign. Since the amount of fluctuation 
is proportional to the temperature, T, and due to the ex- 
clusion principle the average energy of the particle hole 
excitations participating in the rectification effect is pro- 
portional to T as well, pn oc T 2 [jl(|, i.e., it vanishes as 

The second mechanism contributing to p is the quan- 
tum process. In this process, an electron-hole pair cre- 
ated in one of the two layers, tunnels into the second 
layer and is annihilated there. The creation and annihi- 
lation processes occur due to Coulomb interactions with 
the other electrons in the system, which do not tunnel 
between the layers, but which take part in a density fluc- 
tuation that is shared between the two layers as a result 
of the inte rlaye r Coulomb interaction. As we shall see 
in Section IV B below, this process gives rise to correla- 
tions in the momenta of the two layers, in the absence 
of an applied electric field, and hence a mechanism for 
momentum transfer when a field is applied to one of the 
layers. 

To solve Eqs. perturbatively we define: 



V a (r) = V$ ) (r) + V ( « ) (r), 



(12) 



where V^, (r) 



mi i'i ("^(o)^J * S P°t en ti a l (current) 

in the absence of coupling between the layers and 
is proportional to the first order of cou- 



pling between the layers. Lack of coupling to the second 
layer reduces Eq. (||) to 



J$»m=-<r£ i VkV$tf. (13) 

Substitution of Eq. Jl3|)in Eq. (0) (with vanishing right 
hand side for the discussed case) and using Onsager rela- 
tions [see details in Eq. (|24|)] we find the Laplace equation 
AV^q-j = with the boundary conditions (g). The solu- 
tion is straight forward and given by 



(2) 



( 2 )/ 



where 



fr(r)= / d 2 r'p kj (f, 7)^(7) 



(16) 



(17) 



The potential pattern in the second layer depends on the 
location of the pinholes, and the measured resistances 
generally depend on the locations of the voltage contacts. 
Therefore we define integrated potentials by: 

Ui = -l VW(x,y)nidl = - ff ViV^i^r. (18) 
Jds 2 JJs 2 

The last equality in (|l|) follows from Gauss's theorem 

w 



ff V • Gdxdy 



TdS 2 



G ■ ndl with G = wV^ where 



is an arbitrary constant non zero vector. 

If S 2 has a shape of a rectangular of length L 2 and 
width 

W 2 then U x = Jdy [V^ {-L 2 /2, y) - (La/2, y)] and 
U y = Jdx [V^(x,-W 2 /2) - V^(x 1 W 2 /2)]. Thus, the 

x and y components of U are a generalization of the in- 
tegrated transvoltage and Hall transvoltage respectively. 
Using ( |l6| ) we find, in a matrix notation: 



U 



rp^fg\{r) 



s 2 



d 2 rf 2 \r) 



(19) 



s 2 



In order to continue further with the first term we use 



the identity J, 



(2) 



Vj(rijj 2 ') — riVjjy and Gauss's 
7(2) 



(2) 



theorem with Gj — r%Jf - Since p^> does not depend 
on space, we find 



d 2 rpWjW (r) 



S-2 



p^ 1 r 

dS 2 



(n ■ J( 2 )) dl - p& JJ d 2 rrV ■ jgj (r) 



The first term vanishes due to the boundary condition 
(|9~h|) and therefore we find: 



U = P + F, 
where using the Eq. M) we have: 



(20) 



P = -P {2) II d 2 rrJ T (r) = -p^ I/ dhrg'V^ (f) , 
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So 



V, 



(0) 



Jp { xvV, 



V, 



(2) 



(0) 



const (14) an d the vector F is given by 



(i) 
: (o) 



J J (1) - 0- 



d X 



(2) 
(0) 



(2) 
2/(0) 



0. 



(15) 



F = I J "p{r,7)J {1) {7)d 2 rd 2 r'. 



(21) 



(22) 



Treating the coupling between the layers perturbatively, 
the current and the potential in layer 1 arc unaffected, 
to first order, while 



To summarize, we find in this section that the x and y 
components of U are related to the integrated transvolt- 
age and Hall transvoltage and have two contributions. 
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The first term, P describes the leakage contribution to 
the transvoltages. The second term F describes the mo- 
mentum transfer from the first layer to the second due 
to the classical drag effect and due to the quantum con- 
tribution that exist only in the presence of tunneling. 

We have not actually assumed that the magnetic field 
is weak in this section. However, when the Hall angle 
becomes very large, the assumption that the current den- 
sity is uniform in layer 1 may become inappropriate for a 
fixed experimental geometry. In particular, if the length 
L, is not much larger then W, current maybe concen- 
trated near the the edges of the sample, if the Hall angle 
is very large. 



III. DEFINITIONS AND CALCULATION OF THE 
TRANSRESISTANCE AND THE HALL 
TRANSRESISTANCE 

A. F: Momentum transfer due to frictional forces 
and quantum effects 

In the solution of Eqs. (@)-(|) given in Eqs. ©-(H) 
we have assumed that the couplings to the second layer 
(due to frictional forces and/or tunneling) are weak and 
treat them perturbatively. In that approximation p has 
only off diagonal elements (in the layer index) given by 



(/3) 
ij 



(23) 



Notice that in the last equation the layer indices are in 
parentheses to emphasize that here there is no summa- 
tion over these repeated indices. 

For simplicity we assume henceforth that the layers are 
identical and isotropic; extension to non identical layers is 
straight forward. We may then use the following notation 
for different parts of the conductance tensor a 1 * (r, r*): 



a if(r, r) = (TijSap + dij (1 - 5 a p) 



(24a) 



where the intralayer conductivity tensor maybe approxi- 
mated as local and independent of position: 



a a H 
-<j H a 



5(r 



and 



L (r,r 



(24b) 



(24c) 



The conductance a is given by the inverse of the (bare) 
sheet resistance i?g = \/e 2 vD. The Hall conductance 
coefficient <jh is essentially given by Rh/R 2 ^, where 
Rh = H/ne, n is the carrier density, e is the carrier 
charge, and H is the magnetic field. There are correc- 
tions to a and ur due to weak localization and Coulomb 
repulsion in combination with disorder p0| that should 
be included. The tunneling conductance - ^ has, similar 



to <7, corrections due to the interplay between Coulomb 
repulsion and disorder, (see details in Sec. [V). [The cor- 
rection to a and an due to drag and interlayer tunneling 
are small and will be ignored.] 

We have to discuss now the behavior of 5^ (r, f 1 ) . In 
the approximations we use, both frictional forces and 
tunneling amplitudes are small and we obtain to first 
order in the tunneling rate and frictional forces: 



> s 6{? 



+ 0[{<r u /ay,\t\*,\m<r U M]- (25) 



The coefficient a® is related to the drag coefficient pu in 
clean systems (d < I) [O by: 



= <J D S p D /R 2 D , 
ag = (J% = 2a H a D /a = 2R H p D /Rl 



where 



PD 



h C(3)vr 2 1 (T_ 
e 2 16 { K d) 2 {k F d) 2 \~E^ 



(26) 



(27) 



[In Eq. (|26|) we have neglected weak localization 
corrections. If d > I then £[3]/16(fci?G?) 2 — > 
\og[Dn/2dT]/12(kp I) 2 . When correlations between the 
disorder in the layers are included, pn may be multiplied 
by an 0(1) factor 

The term Q in Eq. (^5|) arises due to the interplay 
between the interlayer tunneling and the Coulomb inter- 
action. We will see that the essential contributions to this 
term are from frequencies larger than the temperature, 
i.e., it is of quantum origin, (therefore we use the letter 
Q to denote it). We find in Sec. that Q l jk {r, f) has a 

range cx L T = y/D/T which increases for T — > 0. In the 
expression for F we are interested only in the integral of 
er(r, r 1 ) and the diagonal part of the matrix 



Q« = 



(28) 



was calculated in Ref. |15| (see also Sec. IV Bl) and is 
given (to first order in tunneling) by: 



Q — Qza; — 



e 2 log(rtd) 1 Rq 
~H Kd 2A7rR~' 



(29) 



(Notice that there is a difference of the sample area 
factor divided by L\ between the expression here and 
the expression in Ref. jl5|]. (We introduce it to make 
the temperature dependence clearer.) This expression is 
valid as long as higher order corrections are small and 

In Sec. [VB 2 in Eq. [52| we show that for a weak mag- 
netic field 



Qh = Qxy — —Qyx — 0. 



(30) 
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Using now the definition of F in (|22j) with ([L5|) we find 
to first order in tunneling and frictional forces at weak 
magnetic field 



F x 
F„ 



-JSi n tpD 



J R^QL 2 ^ 
2JR a R H QL 2 T: 



(31) 



where S- ln t — Si S2 is the overlapping region of the re- 
gion of layer 1, Si, and the region of layer 2, Notice 
that in the absence of tunneling F y = 0. The quantum 
corrections oc Q to F x and F y are related by a factor 
2Ru/Rh] this is the analog of the "rule of two" correc- 
tions to the Hall resistance in a single layer samples [E0| . 



B. P: The leakage contribution 

The leakage contribution to transvoltages depends cru- 
cially on the distribution of the pinholes. When the pin- 
holes are distributed uniformly inside a rectangle of size 
a x b centered at the origin we find 



J 



R\ab 



P 



,a/2 






dx 


/ -a/2 




J 


\a 2 



b/2 



-6/2 



dypr [xja + yan / c 2 



12R ± 



b 2 ) a -§y 



(32) 



[Notice that to fulfill the requirement JJ S ^ d 2 rV-J^ = 0, 
we have for a symmetric distribution of pinholes around 
the origin, vij? = const = in Eq. fll4|).] 

The quantum contribution depends on the tempera- 
ture. Thus, to make the analysis complete we should in- 
clude also in the leakage contribution temperature depen- 
dent corrections to R± , a and <Jn arising from intralayer 
electron-electron interactions and weak localization cor- 



rections (see Sec. IV A). Including these corrections we 
find: 



P 



J a* It 



12R ± 



[l + a t t D log £p] 



j(a 2 +b 2 )R H R D F1 . H . 1 11- 

2awi + 2a; n -/3 z ba 



(33) 



where a t = 2a wl + 2ai n -/3zba, a* = 3a in +a wl -/3 z b a and 
<□ = R^e 2 /2ir 2 h. The coefficients a w i, a;„, /3 z b a defined 
in Eqs. ( |37| ) and ( |38| ) are associated with weak localiza- 
tion, intralayer Coulomb repulsion, and Coulomb block- 
ing (Zero bias anomaly) of the tunneling between the 
layers. In zero magnetic field the weak localization cor- 
rections are cut off by the dephasing time, (that is pro- 
portional to 1/T in 2D). In case where H ^> Kc/Dt^ we 
have to cut off the weak localization corrections [ oc a w \ 
in expression 0)] by Th = eHD/hc. 



C. Definition of transresistances 

Now we are at the position to define transresistances. 
As was demonstrated, U x and U y present a generalization 



of transvoltage and Hall transvoltage that are integrated 
on the boundary. In order to define the transresistances 
properly we divide Ui by the current injected from the 
leads to layer 1 and by the appropriate circumference. 
This gives the following definition: 



Rti = 



Ui 



ff s ViVW(x,y)dxdy 



If \m\dl/2 JWi/asJ"*!*/ 2 



(34) 



A rectangular shape of both layers (of size L\ x Wi and 
L2 x W2) reduces these definitions to: 



Rt 



Rt 



U x 



JW1W2 



1 



W 2 I 
and 



1 

L2I 



w 2 /2 



w 2 /2 



dy \v™ (-L2/2, y) - V™ (L 2 /2, y)] (35) 



R 



ty 



l 2 /2 



dx 



L2/2 



JW 1 L 2 

V {2) {x, -W2/2) - V {2) {x, W 2 /2) 



(36) 



(The integration over the edges can be done experimen- 
tally by connecting metallic measuring probes through 
high barriers.) Using the relation U = P + F [see 
Eq. J2C])1 and expressions (|3^) for P and (|3l]) for F, we 
find the expressions for the transresistance (Q) and the 
Hall transresistance (0) cited in the Introduction. 



IV. MICROSCOPIC THEORY OF a, a H , cr D , v%, 

g l , Q l AND Q l H 

In the section we review the dependence of the mi- 
croscopical, geometry-independent parameters a, an, 
cr D , (7^ and g l on the temperature, the interlayer tun- 
neling amplitude, the level of disorder in each layer, 
and the strength of the intralayer and interlayer inter- 
action. We derive the behavior of Q^Ar, f) introduced 



in Eq. (E5J)and derive Q- defined in Eq. fl28|). 



A. Review : g l , a, an and a D 

The effect of disorder and Coulomb interaction in a sin- 
gle 2D layer was studied intensively at the beginning of 
the 1980s p3[|. The interplay between weak disorder, in- 
terference and Coulomb interaction leads to various log- 
arithmic corrections to the sheet conductance a and to 
the Hall conductance a H . It was found that [E3|: 



a = {1/Ra) [1 - i n (a int + a wJ ) log(l/rT)] 



(37) 



where t\j = R.^e 2 /2ir 2 h, and the O(l) dimensionless co- 
efficients, ai n t and a w i, describe the corrections due to 
interaction and weak localization respectively. 
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The parameters g describe tunneling between the two 
layers. The finite conductance inside each layer reduces 
the speed of the charge spreading after the electron tun- 
nels, and leads to a suppression of the tunneling rate 
jl7],|4]]. To first order in the intralayer e-e interaction 

g l ^g l [l-t D {3 zha log(l/TT)], (38) 

where /3 z ba = log(«;d), k = 2"ne 2 v is the inverse Thomas - 
Fermi screening radius. There are no interaction correc- 
tions to a - , however weak localization corrections exist 
and they are twice as larger as the weak localization cor- 
rections to the conductance, 

a H a H [1 - 2t n owi log(VrT)] . (39) 

They tend to be suppressed by a weak magnetic field, 
however in the limit when the magnetic length */ c/2eH 
is larger than the dephasing length ( that is proportional 
to the temperature in 2D) they should be included. 

The expressions for a D and er|^ = 2aH<yo/a are given 
in Eqs. ( |26| ) and (p7f). They were derived in that form in 
Eqs. (35) and (26) of Ref. gl). 

B. The matrix Q\j{r,r) 

By definition the matrix Q\j(f f ,r l ) is given (using a 
linear response formalism) by the correlation function of 

the currents j\ ■ (f f ) with Jj 2 \r'). The Hamiltonian of 
the system is given by [jl5|| 

H = H 1 +H 2 + H int + H T , (40) 

where -Hi(2) is the Hamiltonian of the isolated layer 1 (2), 
including elastic disorder, and 7Jj n t includes interlayer as 
well as intralayer Coulomb interactions. The first three 
terms on the r.h.s. of Eq. ( |40| ) are traditionally involved 
in the description of the drag effect g(]-[l2|. We add 
the term describing pointlikc tunneling processes through 
pinholes 

N 

Ht = E yl E e ii?i ' ( ^4&# + h - c - ' ( 41 ) 

where a(a') and b^) are the annihilation (creation) op- 
erators of electrons in the first and the second layers re- 
spectively, and R l , I = 1 . . . N, are the positions of N 
local pinholes. The coupling energy V 1 is related to the 
tunneling amplitude t by t l = nv\/2LiWiL2W2V l . 

The actual calculations of current-current correlation 
functions are very similar to the calculations of interac- 
tion corrections to the conductivity and the Hall coeffi- 
cient in a single layer. In the description of the calcu- 
lation here step by step we will give references to the 
equivalent steps in case of single layer calculation that 
are described in some details in Ref. Eq]. 



1. Diagonal elements: Q l xx (r,r") 

The largest contributions to Q l xx {r 1 r 1 ) (to second order 
in tunneling amplitude, and first order in interlayer inter- 
action) is coming from the diagrams depicted in Fig. |^, 
(two additional diagrams where the direction of the ar- 
rows is inverted should be added). These diagrams are 
equivalent to the "three diffuson" diagrams contributing 
to the first order intralayer interaction corrections to the 
conductivity of a single layer. [Diagrams 5 (d) and (e) 
of Ref. |2(|.] In the present case, due to the tunnel- 
ing between the layers, the corrections (to second order 
in tunneling amplitude) are due to interlayer interaction 
and there are two diffusons in each layer, i.e., four diffu- 
sons in total, and we expect more singular temperature 
dependence. 



R 1 R 1 




FIG. |2| Two diagram describing the transconductivity matrix 
Qa that are second order in interlayer tunneling (denoted by 
x ) (two additional diagrams with opposite arrows should 
be included). Solid lines with arrows are electrons prop- 
agators. Full circles denote current vector vertexes at 
point r in layer 1 where the external electric field is ap- 
plied and at point r* in layer 2 where the current is mea- 
sured. The crosses denote tunneling through a pinhole 
located at R l . Since impurity (dashed) lines are local 
on scale of the mean free path, the interlayer interaction 
(wavy) line is of range of the order of the distance be- 
tween the layers d, and the relevant frequencies in the 
diffuson (group of dashed lines) are of the order of the 
temperature T, Q xx (f,f^) decays strongly when the dis- 
tance between Ri and f or r* is larger than the thermal 
length L T = \Jd]T . 

These diagrams are calculated using the Matsubara 
frequency formalism |27]] . The analytical structure of the 
frequencies in them is identical to the structure in case 
of single layer interaction corrections to the conductivity 
p6| . In case of a single layer there are, in addition to 
diagrams equivalent to the ones in Fig |2|, "two diffuson" 
diagrams that eventually cancel each other out. [Dia- 
grams 5 (a) (b) and (c) in Ref. ptSj.] This cancellation 
involves a single impurity line [Diagram 5 (c) in Ref. [p6[ ] 
that is not present here, because by assumption there are 
no common impurities (besides the tunneling impurity) 
to the two layers. However, in the present case, due to the 
integration over the different fast momenta in the current 
vertexes and the singular behavior of the four diffuson di- 



8 



agrams, the contribution of two diffuson diagrams is less 
singular and can be neglected. 

After integration over the fast momenta, the formal 
expression for the longitudinal term QL x (f — R l , r' — R l ) 
associated with the diagrams in Fig. g reduces to |L5l : 



Q = E = E J2 J d 2 sd 2 s'Q l xx (s, = 



e 2 1 In(nd) R a 
~%~2Ak nd Rl' 



(48) 



OO 



oj coth 



Q xx {r-R l ^-R 



2T 



F xx {r-R l ,r' -R l -u). (42) 



The function F xx (r — R l , f 1 — R ; uj) to second order in 
the tunneling amplitude is given by: 



FtHr-R 1 ,^ -R l ;uo) = 8 



Rg D 

R l , (L 1 W 1 L 2 W 2 ) 



D(Qx + q x /2)(Q x + k x /2)U 12 (Q, w) x 

Qkq 

(DQ 2 - iuj)- 2 [(D(Q + k) 2 - itu)(D(Q + q) 2 - ioj)}- 1 x 



where 1/R± = J2i V^i- 

The formulas assume that Lt is small compared to 
the system dimensions L and W. When the temper- 
ature T becomes so low that Lt is greater then both 
L and W the transresistance should be flatten out and 
become independent on the temperature. For a system 
where L ^> W, there should also be an intermediate 
temperature regime L > Lt > W, where the system 
becomes quasi-one-dimensional, the integration over the 
momenta in Eq. ( f43| ) are one dimensional and the diver- 
gences are more singular in the temperature. In that case 
integration over the momenta in Eq. (^) [assuming that 
^min > d and therefore the screening properties have a 
2D character] gives: 



(43) 



e iq-(r-R l ) e -ik-{r'-R l ) 



Q = -c- r 



1 Lt 



(49) 



where R l ± = 2irH/e 2 \t l \ 2 is the resistance due to tunneling 

through the pinhole located at ic, and Ui2(Q,oj) is the 
interlayer screened Coulomb interaction, which is given 
in the diffusive case, by pdj |: 



1 



ire 2 q 



(Dq 2 



iuj 



Sint k 2 smhqd V Dq 2 



(44) 



where the divergences at small q are cut off at Dq 2 
u}/(nd). 

Performing the integration over Q we find 



J d 2 sd 2 s'Q l xx {s^) 



1 \n(nd) R a 



h 24tt nd R 



Li 



T t- 



(45) 



where c = 1/ (2tt 7 / 2 ) C (5/2) - 1/13. 

In addition to the exchange diagrams depicted in Fig. || 
diagrams of Hartree type should be included. However, 
unlike the situation in a single layer system p8| , the con- 
tribution of Hartree type diagrams is not singular as the 
contribution of the exchange diagrams of Fig. |^. For- 
mally it happens because the Coulomb line in a Hartree 
diagram does not transfer momenta between the layers. 

Let us examine now what happens for higher order 
terms, i.e., when we have more then a single event of 
interlayer tunneling. For two interlayer tunneling events 
(fourth order in the tunneling amplitude) at point R l 
and R l , we have to consider several diagrams. A typical 
contribution of them is the term: 



On the other hand since the frequencies are of the order of 
the temperature and the momenta Qqk are of the order 
of ' D/uj the function F a and therefore Q l decays at 
range of order of Lt- Physically, the relevant time for the 
quantum phenomena under discussion has to be smaller 
than H/T. At time oc 1/T the electron diffuses for a 
length Lt, hence the quantum correction has a range 
Lt- This allows us to write Q l as 



F^(r~R l y-R l ;r~ * r * 
'Ra D 



r — R l ,f -R l uj) 
2 1 



L 2 W 2 



QUry)^QU(r-7)S l T (\r 



(46) 



where S T (r) is a function of range Lt around R l . We 
normalized it in such a way that / d 2 rS l T {r) = L\. Inte- 
grating the right hand side of Eq. © with respect to r 
and f on all space we find that: 



= L; 



(47) 



Using Eqs. (^H47j, ) we find 



_R± s/LlW^Wi. 
E D(Q x + q x /2)(P x + k x /2)U n (Q,co)x 

PQkq 

(DQ 2 - iuj)- 2 [{D(Q + q) 2 - iuj){DP 2 - iuj)}- 1 x (50) 

(D(P + k) 2 -iujy 1 x 

e iq-{r~R l ) e -ik-{r' e -i(Q-P)-R l e i{Q-P)-R 1 ' _ 

If we assume that the tunneling is at one point than 
integration over the direction of the momentum P causes 
the integral to vanish. When we average on the position 
of the tunneling points at an area axbwe find for a <C Lt 
that the result is less singular than 1/T and for a > Lt 
we recover the results 3(b) and (15) of Ref. [jl5| in which 
Q cx 1/T 2 log T. Other diagrams give similar results. 

Practically, if we are interested in situations when the 
quantum contribution is dominant, we are interested in 
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the case a <C Lt, and we can ignore the second order 
contribution. 




FIG. ^ The most divergent diagrams contributing to the third 
order tunneling (6th order in tunneling amplitude) transcon- 
ductance. Two additional diagrams with opposite arrows 
should be included. 

The third order contributions, involving three tunnel- 
ing events (6th order in tunneling amplitudes), may also 
be computed. The most divergent contributions, de- 
picted in Fig. [| are found by inserting four additional 
crosses in the diagrams of Fig. || in a way that gives the 
maximal number of diffusons. For the case a <C Lt we 
find that their contribution to Q is 



Ro , 1 

log 

nR± Tt 



(51) 



where Q( a ) is the first order contribution to Q given by 
(fl8|). Thus, the first order results are valid as long as 
[Rn/(nR±)] log(l/TV) < 1. In the quasi ID case the 
condition is (Ra/ R±)(Lt / L min ) <C 1. 

By examining carefully the analytical structure of the 
diagrams in Fig. ^, we can try to understand the physics 
leading to the singular quantum contribution (48) to the 
transconductance. 

The ground state of the interacting system may be 
built up from virtual particle-hole (p-h) excitations, rel- 
ative to the Fermi sea of the non-interacting disordered 
system. The e-e interaction, at some instant of time, can 
create two p-h pairs, which propagate forward in time; 
it can annihilate two p-h pairs which originated at an 
earlier time; it can annihilate one p-h pair and create 
another; it can lead to scattering among existing holes 
and/or particles; or it can create or annihilate a single 
pair while scattering an existing particle or hole. The 
most commonly considered correction to the groundstate 
energy of the Fermi liquid, beyond the Hartree-Fock, is 
the RPA correction, which is represented diagrammati- 
cally by closed chains of two or more p-h bubbles. In 
these diagrams, each particle is annihilated by the same 
hole that was created with it originally, and there are no 
scattering processes for existing electrons or holes. 

If one ignores the momentum vertices, the diagram in 
Fig. U is a contribution to the RPA ground state energy 
in which one of the p-h pairs initially created by the 
Coulomb interaction at a time to, tunnels at from one 
layer to the other, before being annihilated at a later time 



t f . The second p-h pair does not tunnel between the two 
layers, but forms part of a polarization cloud, represented 
by the screened interaction propagator, which propagates 
in both layers because of the interlayer Coulomb interac- 
tion. 

Now let us insert momentum vertices in the diagram, 
as shown in Fig 2, so that we can compute the time- 
dependent correlation function between the momenta in 
the two different layers. From this correlation function, 
using the Kubo formula, one may compute the transcon- 
ductance. If the p-h pair is created at time to with to- 
tal momentum q, the polarization propagator will carry 
wavevector —q. Although the system contains impurities, 
so that the electron momentum is not conserved, the di- 
agrams which include averaging over impurity positions 
preserve the total momentum of the p-h pair, and of the 
screened electron propagator. The tunneling processes 
randomize the momenta of the particle and hole, when 
they cross from one layer to the other. Consequently, in 
the absence of Coulomb interactions, there would be no 
correlation in the total momenta of the two layers. Now 
however there is an induced correlation: in order for the 
p-h pair to be annihilated at time tf by the polarization 
propagator created at time to, it must have the same 
total momentum q as the original pair. 

The overall sign of the correlation depends on the sign 
of the interlayer e-e interaction. Indeed, if one changes 
the sign of the interlayer e-e layer, while keeping the fixed 
the interaction between electrons in the same layer, the 
contribution of the diagrams in Fig.^j will simply change 
sign. (This follows from the fact that the bare inter- 
layer interaction appears an odd number of times in the 
screened interaction propagator for this case.) For the 
actual situation, where the interlayer interaction is re- 
pulsive, we find a negative contribution to the transcon- 
ductance. 

When the layers are in the diffusive regime, the parti- 
cle and hole of an excitation stay closer to each other in 
position space, compared to the ballistic case, and stay 
closer to their initial position. This increases the proba- 
bility that they tunnel through the same pinhole, and 
that they eventually annihilate the polarization cloud 
that was originally created together with them. In the 
diagrams, this effect is presented by the dressing of the 
vertexes by the diffuson propagators. 

Our analysis shows that the correlation of the momen- 
tum fluctuations, at T = 0, extends for very long times, 
giving rise to a singular contribution to the transconduc- 
tance. At finite temperatures, when the discontinuity at 
the Fermi energy is rounded, the contribution of the vir- 
tual processes will be cut off at times of order 1/T, and 
the singular contribution is reduced. 
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2. The Hall coefficient Q l xy (r,f*) 

In the absence of a magnetic field and without any e— 
e interactions the electrons' wave front propagating in 
layer 2 (after tunneling from layer 1) is spherically sym- 
metric p5j . This argument is not changed in the presence 
of an external magnetic field, (though the current distri- 
bution has now a component perpendicular to the radial 
direction), and leads us to the conclusion that without 
e—e interaction the tunneling itself between the layers 
does not lead to a finite Hall transconductivity coeffi- 
cient. Generalizing the calculation of the Hall coefficient 
in a single layer [see Figs. 2 (a) and (b) in Ref. ^(|] to 
the case of two layers we find that the diagrams describ- 
ing the Hall transconductivity coefficient, Q l xy = Qrr, 
without e-e interactions, are the ones depicted in Fig.|3|. 
An average over the current direction in the vertexes (full 
dot in Fig. |) proves that this contribution vanishes. The 
insertion of Cooperons [compare to Fig 3. in Ref. p6| ] 
will not change that result. 



OD 



Fig H we find three possible ways to enter the magnetic 
vertex as shown in Fig. [5]. 




FIG. || Diagrams describing the Hall transconductivity that 
are second order in interlayer tunneling (denoted by x) and 
does not include e-e interaction. Dashed vector signifies an 
external magnetic field, solid lines are electron propagators 
and the full dots are current vertexes. An average over the 
current direction at the vertexes shows that this contribution 
vanishes. 

We shall now show that that inclusion of interlayer 
Coulomb interactions does not affect the results that the 
tunneling give no contribution to the Hall transconduc- 
tivity. As in the case without interaction, we insert a 
magnetic vertex in all possible ways to the diagrams con- 
tributing to the transconductivity. Concentrating first in 
the part related to layer 1 of the diagrams depicted in 






FIG. ^ The different ways to insert a magnetic vertex (dashed 
line with an arrow) to the left part of the conductivity cor- 
rections depicted in ^[ These diagrams are similar to those 
of Fig. 6 in jifj. The sum of these diagrams is zero. The 
contribution of "two diffuson" diagrams, (without the diffu- 
son connecting the two sides of the current vertex) vanishes 
as well. 

These diagrams are similar to those of interaction cor- 
rections to the Hall coefficient in a single layer [figure 6 
of Ref. p6|]. It was shown in Ref. |2(| that they cancel 
each other out. 

We showed that there is no noninteracting Hall 
transconductivity and that the most singular terms con- 
tributing to the transconductivity do not give rise to 
any Hall transconductivity as well. Since we look for 
a term different from zero we have to check also that the 
"two diffuson" diagrams does not contribute to the Hall 
transconductivity. But, the insertion of a magnetic field 
vertex to the two diffuson diagrams gives three diagrams 
similar to the one depicted in figure For the "two dif- 
fuson" diagrams the diffuson connecting the two parts of 
the current vertex is missing and the Green functions in 
the two sides of the current vertex has opposite imagi- 
nary part. It can be shown, however, that in this case 
as in the case of the diagrams in Fig. |5| the sum of the 
three vanishes. The calculation is almost identical to the 
calculation in Ref. [^6| and we will not repeat it here. 

This leads us to the conclusion that to first order in 
interaction 

<4 = Q 1 h = Q 1 h = o (52) 

This implies, in turn, that the quantum process gives 
no contribution to the Hall trans conductance. When the 
conductivity matrix is inverted, however, this leads to 
a nonvanishing contribution to the Hall transresistance, 
following the analysis of Sec. [II A. 



V. THE POTENTIAL IN LAYER 2: V (2) (x,y) 

In the previous sections we have discussed the relation 
between integrated potentials measured on the edges of 
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layer 2 and the applied current in layer 1. We have used 
Gauss's theorem and did not need to calculate the po- 
tential in layer 2, V^(x,y) in details. Nevertheless, it 
is instructive to understand, in few simple examples, the 
detailed behavior of V^(x,y). In principle, with ad- 
vanced technology, V ^ (x, y) can be measured directly 
experimentally p9|. We will discuss the two situations 
mentioned in the introduction: the "parallel strip" and 
"cross" geometries with pinhole distribution that is con- 
centrated in the middle of the sample, and identical layers 
with pinholes that are distributed uniformly. In all cases 
we assume that the 2D gases have rectangular shapes of 
sizes L\ x W\ and L 2 x W 2 . 



A. Tunneling through pinholes in the middle of the 
sample 

Now we assume that there are TV pinholes distributed 
inside a rectangular of size a x b at the origin and a,b <C 
Li, L 2 , Wi, W 2 . For simplicity, as before, we assume that 
the frictional forces and the tunneling between the layers 
are weak, i.e., po <C R\j and Rn <C R±- In that case 
we can solve the equations (Q), (@) with the boundary 
condition ([)]) perturbatively in pu/R\j and R\j/Rj_. 

If W 2 < W\ and L 2 < L\ then the classical drag is 
present all over layer 2 and its contribution to the poten- 
tial in layer 2 can be readily found by changing Eq. jl^ ) 
to: 



D vafll 



(53) 



where the diagonal elements of the matrix X a @ are zero 
and the off diagonal are 1. Substitution in Eq. (^) (with 
vanishing right hand side at the discussed case) and using 
(24t) we find the Laplace equation = with the 

boundary condition @ . The solution is straight forward 
and in the limit of <jd *C u given by: 



(i) 



(0) 



J- 



J- 



'H 



T 2 f' "(0) 



JpDX 



(54) 



Notice that the Hall transvoltage vanishes. This happens 
due to the relation @, Q. 

In order to find (perturbatively) the voltage in layer 2 
in the presence of tunneling we substitute Eqs. ( |l2|) and 
( [46| ) in Eq. (||) , keeping terms up to first order in Rjj /R± 
and pd/R± we find: 



J?{r) = J(o)?(f) 



'J 



(55) 



By assumption the corrections to the currents in layer 
1 are small and we neglect them. A subs titution in the 
continuity equation (Q) and a use of ( 24b ) yields the fol- 



(2) 

lowing equation for VA^ , representing the correction to 
the voltage in layer 2 due to tunneling 



>)AV, 



(2) 
(0) 



(0) 



(56) 



The first term in the left hand side of that equation van- 

(2) 

ishes, and we are left with an equation for VA, . The 



(2) 

boundary condition for VA/ should be such that the 

(2) 

residual current due to the potential VA^ vanishes on 
all boundaries, i.e, 

J$ x = aV x V$(±L 2 /2,y) 



J, 



(2) 
(l)t 



-<7 H V y V$(±L 2 /2,y) = Q 
-<t h V x V$(x,±W 2 /2) 
-aV y V$(x,±W 2 /2) = 



(57) 



Using ( p4b| ) for the intralayer conductivities we end up 
with a Poisson equation: 



aAVl 



(2) _ 



S, 



(58) 



where the source term, S, is (— ) the right hand side of 
Eq. ©. 

The Green function of this equation is given by the 
solution of Eq. (^) and boundary condition (|5?]) with a 
source term equal to 5(r — r r ). In the absence of magnetic 
field the boundary conditions can be easily satisfied, by 
introducing a series of image charges. Since the boundary 
conditions are of Neumann type, all the charges have 
the same sign. In the presence of magnetic a field we 
should add to the image charges a source of circulation 
that will cancel the twist of the current field due to the 
magnetic field. The exact forms of the Green function 
G_L 2i w 2 (r, f^) in the cases W 2 <C L 2 and W 2 ^> L 2 are 
given in Eq. (A3). 

For a general source term we find: 

V$(r) = J d 2 r'S(7)G L2 , W2 {ry). 

By assumption, the source term in Eq. j5^ ) is concen- 
trated near the origin and is local, we are interested in 
the voltage at distances much larger than a or 6 and the 
thermal length Lt- In addition, due to the boundary 
condition and conservation of charge J d 2 r' 'S^r") = 0, so 
we can use the a dipole approximation to find: 



V$(r) = P ■ V'G(r), P = J SrrS{r) 
V'G(r) = V'G i2)W2 (f,r%, =0 . 



(59) 



If we assume that the pinholes distribute uniformly 
over a rectangle of size a x 6, the dipole moment associ- 
ated with g l is given by: 



J 



(•a/2 rb/2 

P s = R ^ ab I dx I dy(xx + yy) (R a x + R H y) 



a/2 



6/2 
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p e = 



12R ± 



(a 2 R a x + b 2 R H y) 



(60) 



The dipole moment arising from the quantum corrections 
can be larger than the leakage contribution since at low 
t emp eratures Lj< can be larger than a and than b. Using 
(24t) and the solution (B3) for V we arrive at: 



dxdy x 

]T[xQLV,(^(r)V^ (1) ) x + 
i 

yQLVy (4(r)v x Vo (1) ) y = 

-QJL 2 T (R D x + R H y). 



(61) 



where in the last equality we have used integration by 
parts and the normalization condition J d 2 rS l T (\r\) — 

Lrp . 

Even in the dipole approximation, the exact form of 
V^ 2 \x,y) depends on W2 and L 2 since the information 
about the magnetic field enters through the boundary 
conditions. For the "par allel strip" configuration, i.e., 
L2 ^> W2 using Eq. A4 we finally find, (keeping only 
terms linear in H): 

V 2 (x,y) = Jp D x 
J (_a? 



\Q\L 



2W 2 \12R 

smh{27rx/W 2 ) 2 
cosh(27ra/W 2 ) - cos(27ry/V^ 2 ) ° 



J 
W 2 



a 2 + b 2 



2\Q\L 



12R± 

cosh^x/W^) s\n(-ny /W2) 
cosh(27ra;/Vy2) — cos(27ry/W2) 

sin{2iry/W 2 ) 

2W 2 12R X cosh(27nr/W 2 ) - cos(2iry /W 2 ) 



J 



RoRh 
RuRh 



(62) 



A contour plot of [x, y) with typical parameters is 
depicted in Fig. ^| Since the potentials depend both on 
x and y, the transresistances depend on the locations 
where the potentials are measured. Therefore we have 
defined the integrated transresistances in Eqs. (||) and 
(||). In the "parallel strip" configuration far away from 
the tunneling points, e.g., at x — > ±£2/2, the potential 
practically does not depend on y. Thus, in the "parallel 
strip" configuration while measuring R t the integration 
on the potential is unnecessary. 

In the "cross" configuration where L2 <C W2, we use 
Eq. ( |A5| ) with W2 — > L2, now x and y change their role. 
Far away from the tunneling points, e.g., at y — > ±W2/2 
the potential practically does not depend on x. In that 
"cross" geometry integration on the potential on the 
boundary is not needed while measuring R^ . 




x[\im] 



v'fx,y) [nV] 



FIG. I A contour plot V (2) {x,y) according to Eq. (||) at 
T = 0.5K° . We have used the same sample parameters as in 
the introduction: W2 = 5/irn, mobility /1 = 5 x 10 4 cm 2 /Vs, 
electron density n = 4 x 10 10 cm -2 , R± — 20kQ, a — 
0.1/j,m,b = O.lfim, Rq = 3fcfi, nd = 3, with total current 
I = Q.ly.A and magnetic filed H = 0.05T. 



B. Tunneling through uniformly distributed pinholes 

Now we consider the case where the pinholes are dis- 
tributed uniformly in the sample. In this situation, we 
do not restrict ourselves to first order in the tunneling 
conductance R J . However, we assume that W\ = W% = 
W, L\ = L2 = L. In that case we can approximate 



J dxdy (1/ R±LW ■ ■ •), and substitute the 



expression due to quantum correction in Eq. (|55|) by: 



LW 



(63) 



in that case Eq. (||) should be read: 

J? (f) = [a l3 5 a0 + + o-g) A" Q/3 J VjVg (r) (64) 

Now it is convenient to introduce upper and lower case 
variables: 



(i) 



Ji = Ji 
S = a + (a 1 



J, 



(2) 



f - J (1) - J V 

s = a-{a D +a Q ), 
v = 



7(2) 



v = v^ + v {2 \ 



(65) 



v w - v {2) 



Using now Eqs. (Q), j6^), boundary conditions (^), defi- 
nitions ( |65| ) and the antisymmetry of S and s we arrive 
at the boundary problem: 



at x 



AV = 

±L/2: -S XX V X V — S xy VyV — J (66) 



at y = ± W/2 : -S yx V x V - S yy V y V = 



13 



Av = 

ctt X — — zb L 1 2 ! S x x 

at y = ± W/2 : - -* 
where the 



VxD 



TO V y i; = 



(67) 



characteristic 



momentum q 2 — 2/(R±LWs xx ). Eq. (|67]) is a partial 
differential equation with boundary conditions that are 
not Dirichlet nor Neumann type. A simple straight for- 
ward solution is not available for general values of the 
magnetic field, indicating the complexity of the current 
flow in layer 2 in the presence of a magnetic field and 
uniform tunneling. 

i7) and ( |6^ ) (for weak 
B. We find, then, 



The detailed solution of Eqs. ( 
magnetic field) are given in App. 



1 



Q 

$2 



1 



-y - 



smh(qx) 
coah(qL/2) 



^thc (qL)g(*JL., q L )q W 
Sxx \L W 



(68) 



where thc(x) and g(u,z;a,b) are defined in Eq. (B6) 
and(B7) respectively in App. |B]. 



VI. CONCLUSIONS 

To summarize, we have studied the effect of local pin- 
holes on the drag coefficient, R t , and Hall drag coefficient, 
Rf , in bi-layer dirty 2D systems at weak or vanishing 
magnetic field. We found that there are three contribu- 
tions to R t : the classical drag, the leakage contribution 
and the quantum contribution. The last two exist only 
when there is tunneling through local pinholes between 
the layers. At low temperature the classical drag van- 
ishes. If a, the characteristic size of the region where 
pinholes exist, is smaller then the thermal length, Lt, 
but Lt is small compared to the overall system size, then 
the quantum contribution oc l/T is dominant. 

The measurement of voltage in layer 2 in a typical drag 
experiment does not allow current to flow perpendicular 
its edges. Therefore, the classical drag does not lead to 
any Hall drag coefficient at low temperatures. In con- 
trast, in the presence of pinholes, we find contributions 
to Rf 1 due to the leakage and the quantum contributions. 
As for Rt, if a < Lt the quantum contribution cx l/T is 
dominant at low temperatures. 

The "topographic map" of the actual voltage in layer 
2 is rather complicated. We suggest therefore to study 
Rt in a "parallel strip" geometry (where the the width 
of layer 1 and layer 2 is smaller than their length) and 
to study R^r in a "cross" geometry where the width of 
layer 2 is larger than it length. In these geometries, the 
voltage measurements are made far from the tunneling 
region and the measured potentials do not depend on 
the precise position of the probe along the boundary. 

The drag and Hall drag measurements give direct in- 
formation on the quantum corrections arising from the 
interplay between disorder interaction and tunneling. 



In this work, we have studied the quantum contribu- 
tion up to first order in the measuring current, the tun- 
neling conductance R J , the magnetic field, and the in- 
tcrlayer interaction. We have also estimated the limits 
for the validity of the calculations. 
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APPENDIX A: DERIVATION OF G L ,w(f,r) 

In this appendix we find the Green's function of the 
Poisson equation in 2D: 



aAG 



L,W 



(r, r ) = 8(f— r ) 



(Al) 



with the boundary conditions 

aV x G L , w (±L/2,y;?) + 
a H V y G L w{±L/2,y;7) 
-a H V x G L , w {x,±W/2;7) + 

aV y G LM (x,±W/2-y) 



= 



= 



(A2) 



Without the boundary conditions flA2| ) GL,w(r,r') 



1/27R7 log | 



In the absence of a magnetic field the 



boundary condition can be easily satisfied, by introduc- 
ing a series of image charges. Since the boundary condi- 
tions are of Neumann type all the charges have the same 
sign. In the presence of magnetic field we should add to 
the image charges an image circulation that cancels the 
circulation of the current field due to the magnetic field. 
If we imagine a source charge situated on a black square 
of an infinite chess board then on the black squares we 
have charges of magnitude equal to the source charge and 
on the white squares we should put charges and circula- 
tion sources. Formally, after including the circulation 
sources we find: 

G L ,w{x,y;x' ,y') = 
X - \og[{x-x' n f + {y-y' m f] 



Ana 



n-r-m—e 



l 

Ana 



_m+n— o 



^log [{ X - X ' n f + {y-y' m f 



H 



t a<7 " arctanf^j 



H 



x' n =nL + x'(-l) n , y' m = mW + y'(-l) m . (A3) 

The summation in the first term is over even m + n and 
in the second over odd m + n. 
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In the limiting case when L ^> W we can take only the 
terms with n — 0. The summation of all sources can be 
performed analytically and yields: 

G w {a, b; x, y) = G^^wia, b; x, y) = 



47TCT L 



log 



o<- + a 
4<j<th 



H 



cosh7r 
cosh7r 



x — a 

W 
x — a 

W 



y 



arctan 



'H 



y + b 

tan7r '^^ tanh n 



W 
y + b 

x — a 



2W 
+const 



(A4) 



The constant is arbitrary since the boundary condition 
are of Neumann type. It can be directly checked tha t this 
Green function obeys the boundary conditions (A2) and 
the only singular point inside the strip is at x', y~. Close 
to it we find the usual logarithm which is a solution of 
the Poisson equation with a delta function source term. 
In the "cross" configuration when W ^> L we have: 



G 



(x 



G L (x,y;x',y') = 
-y,V^> x;x' -> -y',y' 



x'). (A5) 



a Neumann type boundary problem, which we can solve 
by Fourier expansion in x/ L. We find 



Vi 

where 



(x,y) = —thc(qL)g(y,^-;qL,qW), (B5) 

s xx q \L W 1 



thc(x) = tanh(a;/2)/0/2) 



(u,z;a,b) 



sinh bz . 
cosh 6/2 * 

,oc J n nf ,L«/,„ 1M s'mha n bz 



v(x,y) 



J 

Sxxq 

s 



sinhf^a;) 



cosh(gi/2) 
y - the ( q L)g(y,^-;qL,qW 

x \-L W 



(B6) 



2Er<cos[ 7 rn("-l)]iS^ (B7) 
and = a 2 /(a 2 + n 2 n 2 ). That gives: 



(B8) 



Expressions (Bl) an d (|B8[ ) are the solution of the bound- 
ary value problems (|66|) and ( |67|) in the limit of a weak 
external magnetic field. 



APPENDIX B: SOLUTION OF EQS. (|66|) AND 
@ FOR WEAK MAGNETIC FIELD 



The solution of Eq. | 
wardly as 



can be written straightfor- 



V(x,y) = 



J S x 



det 



J $xy 

det 



(Bl) 



where det = det S = S xx + S xy . This resembles the 
solution of the usual single layer Hall problem. 
At zero magnetic field the solution of Eq. (|67[) is: 



vo(x,y) 



J sinh(gx) 



s a;2; gcosh(gL/2)' 
For weak magnetic fields, expanding v around vq, 



v 



^xy 



(B2) 



(B3) 



and substituting in Eq. (67) we find 



A«i = q 2 vi 
at x = ±L/2 : V ' x v Y = 
at y = ±W/2 : V yVl = V x v , 



(B4) 



where we have kept only terms linear in the magnetic 
field H . (Notice that the off diagonal elements of the 
matrices S and s are proportional to H.) Thus in the 
case of a weak magnetic field the problem is reduced to 
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